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ABSTRACT 

Indirect dark matter searches with ground-based gamma-ray observatories 
provide an alternative for identifying the particle nature of dark matter that is 
complementary to that of direct search or accelerator production experiments. 
We present the results of observations of the dwarf spheroidal galaxies Draco, 
Ursa Minor, Bootes 1, and Willman 1 conducted by the Very Energetic Radi- 
ation Imaging Telescope Array System (VERITAS). These galaxies are nearby 
dark matter dominated objects located at a typical distance of several tens of 
kiloparsecs for which there are good measurements of the dark matter density 
profile from stellar velocity measurements. Since the conventional astrophysical 
background of very high energy gamma rays from these objects appears to be 
negligible, they are good targets to search for the secondary gamma-ray pho- 
tons produced by interacting or decaying dark matter particles. No significant 
gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a 
typical exposure of ~ 20 hours. The 95% confidence upper limits on the integral 
gamma-ray fiux are in the range 0.4 — 2.2 x 10~^^ photons cm~^s~^. We interpret 
this limiting fiux in the context of pair annihilation of weakly interacting massive 
particles (WIMPs) and derive constraints on the thermally averaged product of 
the total self-annihilation cross section and the relative velocity of the WIMPs 
{{(jv) < 10^^'^ cm'^ s^^ for > 300 GeV/c^). This limit is obtained under con- 
servative assumptions regarding the dark matter distribution in dwarf galaxies 
and is approximately three orders of magnitude above the generic theoretical pre- 
diction for WIMPs in the minimal supersymmetric standard model framework. 
However significant uncertainty exists in the dark matter distribution as well as 
the neutralino cross sections which under favorable assumptions could further 
lower this limit. 

Subject headings: gamma rays: observations — dark matter — galaxies: dwarf 

1. Introduction 

The existence of astrophysical non-baryonic dark matter (DM) has been established by 
its gravitational effects on a wide range of spatial scales. Perhaps the most compelling evi- 
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dence for the existence of weakly interacting particle dark matter comes from observations 
of colliding galaxy clusters in which the baryonic matter in the form of X-ray emitting gas is 
sepa rated from the source of the gravitatio nal potential detected through gravitational lens- 
ing fjClowe et a/.ll2006l : iBradac et a/.ll2008l ). However, despite the well-established presence 
of DM in the universe, its particle nature is unknown. 

The quest to understand the nature of DM draws upon research in cosmolo gy, particle 
physics, and astroparticle physics with direct and indirect detection experiments (IBergstrom 
2OOOI : iBertone et a/.ll2005[ ). In this paper we focus on the indirect search for very high energy 
(VHE, Energy > 100 GeV) gamma rays resulting from the interaction or decay of DM 
particles in astrophysical objects in which the gravitational potential is dominated by DM. 



Among many theoretical candidates for the DM particle ( iTaoso et a/.ll2008l ). a weakly 
interacting massive particle (WIMP) is among the best motivated. A thermal relic of the 
early universe with an interaction cross section on the weak scale will naturally produce 
the present-day D M density if the particle has a weak-scale mass (ILee fc Weinbergl 119771 : 
D^^Tdlmiiji^DMh'^ = 0.1099 ± 0.0062 (WMAP only), n^Mh'^ = 0.1131 ± 0.0034 
(WMAP -|- Baryon Acoustic Oscillations + Type la Supernovae), where floM is the ratio of 
dark matter density to the critical density for a fiat universe and h is a dimensionless quanti ty 
defined as the Hubble constant. Ho, normalized to 100km s~^Mpc~^ f lKomatsu et a/.ll2009l )). 
Several candidates for WIMPs are predicted in extensions to the standard model of particle 
physics, for example, the neutralino from supersvmmetr v (lEUis et a/.lll984l) and the Kaluza- 



Klein particle in theories of universal extra dimensions ( IServant fc Taitll2003t iBertone et al 



2OO3I ). Both neutralinos and Kaluza-Klein particles are predicted to have a mass in the range 
of a few tens of GeV/c^ to possibly a few TeV/c^. 

The self-annihilation of WIMPs produces a unique spectral signature of secondary 
gamma rays which is expected to significantly deviate from the standard power-law be- 
havior observed in most conventional astrophysical sources of VHE gamma rays and would 
have a cutoff at the WIMP mass. In addition, it could exhibit a monoenergetic line at the 
WIMP mass or a considerable enhancement of gar nma-ray photons at the endpoint of the 
spectrum due to the internal bremsstrahlung effect ( iBringmann et a/.ll2008[ ). Observation of 
these spectral signatures combined with the spatial distribution of the gamma-ray fiux from 
an astrophysical source is a unique capability of indirect DM searches utilizing gamma rays. 

Nearby astrophysical objects with the highest dark matter density are natural candi- 
dates for indirect DM searches. While the Galactic Center is likely to be the brightest 



source of annihilation radiation (e.g. see IBergstrom et a/.lll998l ). VHE gamma-ray measure- 



ments reveal a bright gamma-ray sou rce at the center which constitutes a large astrophysical 
background (lAharonian et a/.ll2006a[ ). Other possible bright sources are expected to be the 
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cores of nearby large galaxies such as M31 or halos around galactic intermediate mass black 
holes, should they exist, where adiabatic compression of dark matter halos could result in 
a large enha ncement in the annihilat i on signal, in some ca ses already exceeding experimen- 
tal bounds (IBringmann et a/.ll2009al : iBertone et a/.ll2009l ). However, in these regions, the 
DM density profiles are poorly constrained and, in the case of nearby galaxies, conven- 
tional astrophysical VHE sources can generate backgrounds for DM annihilation searches. 
In contrast, the satellite dwarf spheroidal galaxies (dSphs) of the Milky Way are attractive 
targets for indirect dark matter searches due to their proximity (20-100 kpc) and relatively 
well- constrained DM profiles derived from stellar kinematics. They, in fact, may be th e 
brightest sources for annihilation radiation after the Galactic Center (iBuUock et a/.ll2009l ). 
The general lack of active or even recent star formation in most dSphs implies that there is 
little background from conve ntional astrophysical VHE processes as has been o bserved in the 
Milky Way Galactic Center JXosack et a/.ll2004i kharonian et a/.ll2006bl.l2009[). The growing 



class of nearby dSphs discovered by recent all-sky surveys (lYork et al 



200 



Belokurov et al 



20071 ) increases the probability of finding an object for which the halo density is sufficient to 
yield a detectable gamma-ray signal. 

In this paper, we report on an indirect DM search for gamma rays from four dSphs: 
Draco, Ursa Minor, Bootes 1 and Willman 1, carried out using the Very Energetic Radiation 
Imaging Telescope Array System (VERITAS). After a brief summary of the properties of the 
observational targets and previous VHE observations in Section 2, we describe the VERITAS 
instrument, the data set, and the analysis techniques in Section 3. Sections 4 and 5 are 
devoted to the discussion of the results and their interpretation in terms of constraints on 
the WIMP parameter space. We conclude in Section 6 with a discussion of the opportunities 
for indirect DM detection by future ground-based gamma-ray instrumentation. 



2. Observational Targets 



Three of the dSphs forming the subject of this paper, Draco, Ursa Minor, and Willman 
1, have been identified as the objects wi thin the dSph class with potentially the highest 
gamma-ray self-annihilation fiux, e.g. see IStrigari et al\ (120071 . 120081 ). The modeling of the 
DM distribution of these galaxies usually is based on stella r kinematics assumin g a spherically 
symmetric stellar population and an NEW profile for DM ( iNavarro et a/.lll997l ) characterized 
by two parameters: the scale radius and scale density ps, 



p{r) = Ps 



r 



(1) 
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Table 1: Properties of the four dSphs. Preferr ed values for DM ha lo p arameters, p., and 



r..;, whic h are defined in the text are taken from IStrigari et all (120071 ) and iBringmann et al. 
(l2009bf). Val u es for Ly, the visual luminosity, and Th, the half-light radius are taken from 



Walker et al\ ( 120091 ) . Rd is heliocentric distance of the dSph. The calculation of the dimen- 
sionless line of sight integral, J, which is normalized to the critical density squared times the 
Hubble radius (3.832 x lO^^GeV^cm"^), is explained in Section |5l The J value for Bootes 
was calculated by G.D. Martinez and J.S. Bullock. As explained in the text, the elongation 
of Bootes and the relative lack of stellar kinematic data lead to large uncertainties for or 
Ps and no values are provided in this case. 



Quantity 


Draco 


Ursa Minor 


Bootes 1 


Willman 1 


a [J2000.0] 


17^20'"12.4" 


15^ 09'^11.3" 


14^00^06" 


10M9'"22.3" 


S [J2000.0] 


57°54'55" 


67°12'52" 


14°30'00" 


5r03'03" 




(2.7 ±0.4) X 10^ 


(2.0 ±0.9) X 10^ 


(3.0 ±0.6) X 10^ 


(1.0 ±0.7) X 10^ 


r/i [pc] 


221 ± 16 


150 ± 18 


242 ± 21 


25 ±6 


Rd [kpc] 


80 


66 


62 


38 


Ps [Me/kpc^] 


4.5 X 10^ 


4.5 X 10^ 




4 X 10^ 


Ts [kpc] 


0.79 


0.79 




0.18 




4 


7 


3 


22 
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The p rop er ties of these ga 



axies including constraints on and ps as found in IStrigari et al 



(120071 ) and IStrigari et al\ (120081 ) are summarized in Table [H 



The Draco dSph is one of the most frequeritly studied objects for ind i rect DM detection 



2007 



(iBaltz et a/.ll2000l:lTvlei 



Strigari et a/.l 12007 . 



2002 



2008 



Evans et a/j2004: Colafrancesco et a/.ll2007: Sanchez-Conde et al. 



Bringmann et g/J l2009br) . It has an approximately spheri 



cally symmetric stell ar distribution (Irwin fc Hatzidimitrioulll995[ l with total luminosity of 
the ord er of 10^ L^r^ (jPiatek et a/.l 120021). The large spectroscopic dat a set available for this 



object (IWilkinson et al\ l2004l : iMunoz et al\ l2005l : IWalker et al\ l2007l ) tightly constrains its 
DM distribution profile. Draco is consistent with an old low-metallicity ([ Fe/H] = -1.8 ± 



0.2) s tellar population with no significant star formation over the last 2 Gyrs (lAparicio et al 
20011 ). Draco p revio usly has been observed at V HE energies by th e STACEE observatory 



/UUll). JJraco previo usly nas been observed at \ niij energies by tn e DiAUijjiij observatory 
jDriscoU et adbooslj. the Whipple 10m telescope JWood et a/.ll2008h . and the MAGIC tele- 
scope jAlbert et a/.ll2008h . 



The Ursa Minor dSph has a distance and inferred DM content similar to those of 
Draco. There is no ev idence of young or intermediate age stellar populations in Ursa Minor 
( IShetrone et a/.ll200ll ). Photometric studies of this objec t have found evidence for significant 
struc tures in the stellar distribution in th e central 10^ (jBellazzini et al. 2002; Kleyna et al. 



20031 ) and an extratidal stellar population (jPalma et a/.ll2003l ). These unusual morphological 
characteristics could be evidence of possible tidal interaction with the Milky Way, velocity 
projection effects along the line of sight, or the presence of fluctuations in the DM induced 
gravitational potential ( IKleyna et a/.l 120031 ). In fact, such confusing factors are present in 
most dSph gal axies. Ursa Minor was previously studied at VHE energies by the Whipple 
10m telescope dWood et a/.1l2008h . 



The recently discovered dSph Bootes 1 (IBelokurov et a/.ll2006l ) shows evidence for elon- 
gation of the stellar profile. N-Body simulations can not reproduce the observed velocity 
dispersion without a dominant contribution from DM. In addition, modeling of the tidal in- 
teraction effects between Bootes 1 and the Milky Way do not provide an adequate explanation 
for the elongation of this system suggesting a n on-spherically symmetric distribution of DM 
in the Bootes progenitor (jPellhauer et q/.II2008| ). Given the stellar kinematical data is based 
on about 30 stars, the scale radius and density of the NFW profile have large uncertainty 
as well as signficant degeneracy. Thus, values for and ps are unavailable in the litera- 
ture. The modeling of Bootes 1 was done by G.D. Martinez and J.S. Bulloc k (2009, private 



comm u nication) for a ran ge of NFW fits. The methodology is described in iMartinez et al. 



(120091 ): lAbdo et al\ (120101 ). They produce a probability density function (pdf) for J, the 
astrophysical contribution to the flux (see Section [5]), which is approximately Gaussian in 
log(J). The value given in Table [1] represents J at the peak of the pdf which is approx- 



- 8 - 



imately the mean of the distribution. The estimates of the age of the stellar population 
and metallicity sugg est similarity with the old and metal-poor ([Fe/H] ^ -2.5 - - 2.1) stellar 



distribution of M92 flBelokurov et a/.ll2006l : iMunoz et a/.ll2006l : Martin et a/.l 120071 ). To date, 
no other VHE gamma-ray observations have been reported for this object. 

Together with Bootes 1, Willman 1 belongs to the ne w class of low surfac e brightness 
dSphs recently discovered by the Sloan Digital Sky Survey (IWillman et a/.ll2005l ). Willman 1 
is one of the smallest {r^ ~ 25 pc) and least luminous (Ly ~ 10^ Lq) dSphs known. Its half- 
light radius and absolute magnitude s uggest that it may be a n intermediate object between 
dwarf galaxies and globular clusters ( iBelokurov et a/.ll2007| ). Due to the small kinematic 
sample of stars available for this object, the constraints on the DM halo parametrization 
are poor ( IStrigari et a/.ll2008l ). Latest estimates of the metallicity suggest a low value of 
[Fe/H] which is con sistent with the observed trend of decreasing metallicity for fainter dSphs 
( ISiegel et al\ |2008| ). The MAGIC collabor ation has recently reported the results of the 
observations of Willman 1 (lAliu et a/.ll2009l ). 



3. Data and Analysis 
3.1. The VERITAS Observatory 



The VERITAS observatory is an array of four 12m imaging atmospheric Cherenkov 
telescopes (lACTs) located at the Fred Lawrence Whipple Obser vatory (31°57^N 111 °37^W) 
in southern Arizona at an altitude of 1.27 km above sea level (IWeekes et a/.ll2002l ). The 
observatory is sensitive over an energy range of 150 GeV to 30 TeV with an energy resolution 
of 10-20% and an angular resolution (68% containment) of < 0.14° per event. For the 
measurements reported here, VERITAS had a point source sensitivity capable of detecting 
gamma rays with a flux of 5% (1%) of the Crab Nebula flux above 300 GeV at five standard 
deviations in < 2.5 (< 50) hours at 20° zenith angle. During summer, 2009 subsequent 
to the four dSph observations, the array configuration was changed, improving the point 
so urce sensitivity. Fur ther technical description of the VERITAS observatory can be found 
in 



Acciari et al.\ (120081) 



3.2. Data 



Observations of the Draco, Ursa Minor, Bootes 1, and Willman 1 dSphs wer e performed 
durin g 2007-2009 (see Table |2]). Observations were taken in "wobble" mode ( iBerge et al. 
20071 ) with the source offset by 0.5° from the center of the field of view in order to obtain 
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source and background measurement within the same observation. The direction of the offset 
was alternated between north, south, east, and west to minimize systematic errors. Reflected 
background regions are defined within the field of view at the same radius with respect to the 
camera center as that of the targeted dwarf galaxy. Observations were made with varying 
atmospheric conditions during moonless periods of the night. Data were quality selected 
for analysis based on the stability of the cosmic-ray trigger rate and the rms temperature 
fluctuations observed by an FIR camera viewing the sky in the vicinty of the observed target 
(< 0.3°C). The total exposure on each source is given in Tabled 



3.3. Analysis 



Data reduction follows the methods described in lAcciari et all (120081 ). A brief outline of 
the analysis flow follows. Images recorded by each of the VERI TAS telescop es are character- 
ized by a second moment analysis giving the Hillas parameters (jHillaslll985l ). A stereoscopic 
analysis of the image parameters is used to reconstruct the gamma ray arrival direction 
and shower core position. The background of cosmic rays is reduced by a factor of > 10^ 
utilizing cuts on the reconstructed arrival direction (6*^ < 0.013 deg^) and the image shape 
parameters, mean scaled width and length (0.05 < msw < 1.16 and 0.05 < msl < 1.36). 
The image distance from the center of the camera is required to be less than 1.43° to avoid 
truncation effects at the edge of the 3.5° field of view. The integrated charge recorded in 
at least two telescopes is further required to be > 75 photoelectrons (400 digital counts) 
which effectively sets the energy threshold of the analysis to be above ~200 GeV depending 
on the zenith angle. The energy threshold quoted in our analyses is taken to be the energy 
at which the differential detection rate of gamma rays from the Crab Nebula peaks. The 
cuts applied in this analysis were optimized to maximize significance of the detection for 
a hypothetical source with a power-law spectrum (dF/dE = 3.2x 10~^^ (E/TeV)~^'^ cm~^ 
s~^ TeV~^) corresponding to 3% of the Crab Nebula flux. Two independent data analysis 
packages were used to analyze the data and yielded consistent results. 



Table 2: Summary of observation periods and exposures of dSphs by VERITAS. 



Source Period Exposure [hr] Zenith Angle 



Draco 2007 Apr-May 18.38 26 - 51 

Ursa Minor 2007 Feb-May 18.91 35 - 46 

Bootes 1 2009 Apr-May 14.31 17 - 29 

Willman 1 2007 Dec-2008 Feb 13.68 19 - 28 
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The significance of tlie detection was calculated by comparing the counts in the source 
region to the expected background counts. The background in the source region is estimated 
using the reflected region model. In this model circular background regions, here of angular 
radius 0.115°, are defined with an offset from the camera center equal to that of the putative 
source. Eleven background regions can be accommodated within the VERITAS field of view. 
The absence of bright stars within any of the four dSph pointings allows all eleven regions 
to be used in the background count estimati on. The significance of any signal is computed 
using the Li and Ma method (ILi fc Malll983l eqn. 17). 



4. Results 



Table [3] summarizes the results for each of the four dSphs. The effective energy threshold 
for each of the targets changes primarily due to the average zenith angle of observations. The 
table shows the average effective collecting area for gamma rays as calculated from a sample 
of simulated gamma-ray showers. No significant excesses of counts above background were 
detected from these observations. The 95% confidence level upper limits on the gamma-ray 
integral flux were calculated using the bounded profile likelihood ratio statistic developed by 
Rolke et all tood ). 



As we have noted, Draco, Ursa Minor, and Willman I have been observed by other 
lACTs and we briefl y compare our flux limits to the other observations. For Draco, STACEE 



( DriscoU et a/.lboosi ) fin ds a spectral limit o f less than 1.6 x 10 ( 200Gev ) cm ^s ^GeV 



s 

-11 



— 2 —1 

cm s 



The MAGIC flux limit fl Albert et fli.1120081 ) from observation of Draco is 1.1 x 10 
above a threshold of 140 GeV. MAGIC also has set flux limits for Willman 1 in the range 
5.7 — 9.9 X 10~^^cm~^s~^ above 100 GeV based on several benchmark models (lAliu et al. 
20091 ). compared to our limit of 1.17x 10~^^cm~^s~^ above a threshold of 320 GeV. The limits 
for the VERITAS observations of Draco and Ursa Minor are an improveme nt of about a facto r 
of 40 over the earlier observations of the group on the Whipple 10m lACT ( jWood et a/.ll2008h . 



Figure [T] shows the upper limits on the differential spectral energy density (E^ d0/dE) 
as a function of energy. The upper limits were derived with four equidistant log energy bins 
per decade requiring 95% C.L. in each bin. 
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5. Limits on WIMP Parameter Space 

The differential ffux of gamma rays from WIMP self-anniliilation is given by 



d4>{An) (av) 



dE Svrm^ 



dN{E,m^ 
dE 



[ dn [ p^{x,n)dx, (2) 

J An J 



wfiere (erf) is tlie tliermally averaged product of the total self-annihilation cross section and 
the velocity of the WIMP, is the WIMP mass, dN{E, m^)/dE is the differential gamma- 
ray yield per annihilation, AQ is the observed solid angle around the dwarf galaxy center, 
p is the DM mass density, and A is the line-of-sight distance to the differential integration 
volume. The astrophysical contribution to the flux can be expressed by the dimensionless 
factor J 

J{An) = (^—] [ dfl [ p\X,Q)dX, (3) 



pIRh / J An 

which has been normalized to the product of the square of th e critical den s ity, P r = 9.74 x 
Xg-30g cm~^ and the Hubble radius, Rh = 4.16 Gpc following IWood et all (120081 ). 



Based on equation [21 the upper limits on the gamma-ray rate, i?^(95% C.L.), constrain 
the WIMP parameter space (m^, (erf)) according to 

R^{95% C.L.) ^ J f {av) 



hr"^ 1.09 X 104 V3 X 10-26cm3s-i 

/•°° A{E) /SOO GeV/c^y EdN/dE{E,m^)dE 
^ Jq 5 X 108cm2 V m'^ ) iF^ 

where A{E) is the energy- dependent gamma-ray collecting area. The expression has been 
cast as a product of dimensionless factors with the variables normalized to representative 
quantities, e.g. the cross section times velocity is normalized to 3 x 10~26cni3s~^ which is a 
rough generic prediction for (av) for a WIMP thermal relic in the absence of coannihilations 
for > 100 GeV/c^ (c.f. figure [2]). The main contribution to the integral comes from 
the energy range in the vicinity of the energy threshold {E ^ 300 GeV for observations 
in this paper) where A{E) changes rapidly. For VERITAS the effective area at 300 GeV 
is ~ 6 X 10*^ cm^. For a representative MSSM model, EdN/dE at 300 GeV is a function 
of neutralino mass, m^, and it changes in the range 10"^ — 10~^ for from 300 GeV/c^ 
to a few TeV/c^. Although EdN/dE is a rapid function of m^, this dependence is nearly 
compensated by the (300 GeV/c^/m^)^ prefactor. The product of these two contributions 
and, consequently, the overall integral value, is weakly dependent on the neutralino mass 
within the indicated range and is on the order of 1. It is evident from the inequality (Equation 
H]) that for a typical upper limit on the detection rate of 1 gamma ray per hour, significantly 
constraining upper limits on (av) could be established if J is on the order of 10^. 
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Because the factor, J, is proportional to the DM density squared, it is subject to con- 
siderable uncertainty in its experimental determination. For example, the mass of a DM 
halo is determined by the interaction of a galaxy with its neighbors and is concentrated 
in the outer regions of the galaxy. Unlike the DM halo mass, the neutralino annihilation 
flux is determined by the inner regions of the galaxy where the density is the highest. For 
these regions the stellar kinematic data do not fully constrain the DM density profile due 
to limited st atistics. Various parametrizations of the DM mass density profile have been 
put forward (INavarro et a/.lll997l : iKazantzidis et a/.ll2004l : Ide Blok et a/.ll200ll : lBurkertlll995l ) 
based on empirical fits and studies of simulated Cold D ark Matter (CDM) halos. We adopt 
the assumption of the NFW profile (INavarro et a/.lll997l ) given in Equation [1] which describes 
a smooth distribution of DM with a single spatial scale factor r^. The astrophysical factor, 
J, is then given by 



27rp: 



pIRh 



COs(0.115 ) J \mi 



r(A) 



1 + 



r(A) 



n -4 



dX d{cos9), 



(5) 



where the lower integration bound of 0.115° corresponds to the size of the signal integration 
region. The galactocentric distance, r(A), is determined by 



r(A) 



-^^ + R-dSph - '^^Rdsph cos 6*, 



(6) 



where A is the line of sight distance and RdSph is the distance of the dwarf galaxy from the 
Earth. 

Although the integration limits, A^m and Xmax, are determined by the tida l radius of 
the dSph {rt = 7 kpc was used for these calculations) (ISanchez-Conde et a/.ll2007l ). the main 
contribution to J{AQ) comes from the regions r < ^ and therefore the choice of 
negligibly affects the J value. The main uncertainty for J computation is due to the choice of 
Pk and r.s. For Draco and Ursa Minor, ps and are taken as the r nidpoints of the range from 
Strigari et al\ (120071 ) . For Willman 1 , and are adopted from iBringmann et all (j2009bl ) . 
The J value Bootes 1 was calculated by Martinez and Bullock as discussed in section |21 The 
summary of the J values calculated for each object is given in Table [TJ 

An estimated value of J of order 10 is representative for all observed dSphs, which is 
three orders of magnitude below the value needed to constrain generic WIMP models with 



> 100 GeV/c^. Figure [2] shows the exclusion region in the 



m 



X' 



(cTf)) parameter space 



due to the observations reported in this paper. MSSM models shown in the figure were 
produced with a random sca n of the 7-parameter phase s pace defined in the DarkSUSY 
package (iGondolo et a/.l 120041 ) with the additional WMAP ( ISpergel et a/.l 120071 1 constraint 
on the cosmological DM energy density. 
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Several astrophysical factors can increase the value of J as compared to the conservative 
estimates given in Tabled] First, the inner asymptotic behavior of the DM density may be 
steeper than oc predicted by the NFW profile due to unaccounted physi cal processes 
at sm all spatial scales. The extreme assumption would be the Moore profile ( iMoore et al 



19991 ) oc r~^'^ asymptotically which generates a logarithmically divergent self-annihilation 
flux indicating that another physical process, for example self-annihilation, would limit the 
DM density in the central regions of the galaxy. A second factor that would increase the 
value of J is deviations of the DM distribution from a smooth ave rage profile (substruc- 



tures). CDM N-body simu lations predict substructures in DM halos (ISilk fc Stebbinslll993 



Diemand et a/.ll2005l 120071 ) and the effects on the DM self-annihilation have been studied in 
these simulations. In general any regions of DM overdensity will enhance the self-annihilation 
flux; the cumulative effect of these enhancements is usually referred to as the boost factor. 
Strigari et al\ (120071 ) find a maximum boost factor of order 10^ while a more detailed cal- 



culation that accounts for the particle properties of the neutralino duri ng formation of DM 
halos suggests boost factors of order 10 and below ( jMartinez et a/.ll2009l ). Thus, present gen- 
eration lACTs could be as close as two orders of magnitude in sensitivity from constraining 
generic MSSM models. 

Two effects related to the properties of the WIMP particle may improve the chances 
of the detection of neutralino self-annihilation by ground-based gamma-ray observatories. 
Internal bremsstrahlung gamma rays produced in neutralino self-annihilation recently calcu- 
lated by lBringmann et al\ (|2008[ ) can significantly enhance dN/dE at the energies comparable 



to niy^ for some MSSM models due to the absence of the helicity suppression factor. Effec- 
tively this increases the value of the integral in Equation HI especially for the higher mass 
neutralino models. In addition, the {av) for self-annihilation at the present cosmological 
time may be considerably larger than at the time of WIMP decoupling due to a velocity- 
dependent term in the cross-section and the reduct ion of the kinetic energy of the WIMP 
due t o the cosmological expansion of the universe (IRobertson fc Zentnerl 120091 : iPieri et al. 
2009[ ). 



6. Conclusions 

We have carried out a search for VHE gamma rays from four dSphs: Draco, Ursa 
Minor, Bootes 1, and Willman 1, as part of an indirect DM search program at the VERITAS 
observatory. The dSphs were selected for proximity to Earth and for favorable estimates of 
the J factor based on stellar kinematics data. No significant gamma-ray excess was observed 
from the four dSphs, and the derived upper limits on the gamma-ray fiux constrain the (erf) 
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for neutralino pair annihilation as a function of neutralino mass to be < 10~^^ cm^ s~^ for 
~ 300 GeV/c^. The obtained {av) hmits are three orders of magnitude above generic 
predictions for MSSM models assuming an NFW DM density profile, no boost factor, and 
no additional particle-related gamma-ray flux enhancement factors. Should the neglected 
effects be included, the constraints on (av) in the most optimistic regime could be pushed 
to < 10-25 cm^ s-\ 

To begin confronting the predictions of generic MSSM models through observation of 
presently known dSphs, future ground-based observatories will need a sensitivity at least 
an order of magnitude better than present-day instruments. The list of dSphs favorable for 
observations of DM self-annihilation has grown over the last years by a factor of roughly two, 
and it is anticipated that newly discovered dSphs may offer a larger factor J. The ongoing 
sky survey conducted by the Fermi Gamma-ray Space Telescope (FGST) may also identify 
nearby higher DM density substructures within the MW galaxy which could be followed up 
by the lACT observatories. Typical current exposures accumulated on dSphs by lACTs are 
of order 20 hours, and ongoing observing programs could feasibly increase the depth of these 
observations by a factor of 10 (a sensitivity increase of ~3). Improvements in background 
rejection are anticipated to increase sensitivity by an additional 20-50%. The soon-to-be- 
operational upgrades, MAGIC-II and HESS-II, as well as a planned VERITAS upgrade will 
reduce the energy threshold and consequently increase the dN/ dE contribution by a factor 
as large as 10 thus providing an additional sensitivity improvement. With all these factors 
combined, the {av) limits for > 300 GeV/c^ will begin to rule out the most favorable 
MSSM models assuming a moderate boost factor. Next generation lACT arrays now being 
planned such as the Advanced Gamma-ray Imaging System (AGIS) and the Cherenkov 
Telescope Array (CTA) will provide an order of magnitude increase in sensitivity and lower 
the energy threshold by factor of ~2 as compared to VERITAS. These instruments will be 
able to probe the bulk of the parameter space for generic MSSM models with ~ 300 
GeV/c^ without strong assumptions regarding potential flux enhancement factors. 

This research is supported by grants from the US National Science Foundation, the US 
Department of Energy, and the Smithsonian Institution, by NSERC in Canada, by Science 
Foundation Ireland, and by STFC in the UK. We acknowledge the excellent work of the 
technical support staff at the FLWO and the collaborating institutions in the construction 
and operation of the instrument. V.V.V. acknowledges the support of the U.S. National 
Science Foundation under CAREER program (Grant No. 0422093). 



^http://www. agis-observatory.org/ 
^http: / / www.cta-observatory.org/ 
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Table 3: Results of observations of dSphs by VERITAS. 



Quantity 


Draco 


Ursa Minor 


Bootes 1 


Willman 1 


Exposure [s] 


66185 


68080 


51532 


49255 


On Source [counts] 


305 


250 


429 


326 


Total Background [counts] 


3667 


3084 


4405 


3602 


Number of Background Regions 


11 


11 


11 


11 


Significance'^ 


-1.51 


-1.77 


1.35 


-0.08 


95% C.L. [counts]*^ 


18.8 


15.6 


72.0 


36.7 


Average Effective Area [cm^] 


5.84 X 10^ 


5.71 X 10^ 


6.37 X 10^ 


6.37 X 10^ 


Energy Threshold [GeV]*^ 


340 


380 


300 


320 


Flux Limit 95% C.L. [cm'^s"^] 


0.49 X 10~i2 


0.40 X 10-12 


2.19 X 10-12 


1.17 X 10-12 



" Li and Ma method (|Li fc Malll983l) . 
''Rolke method |Rolke et alWzm^ . 
'^Definition given in text. 
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Fig. 1. — 95% CL. upper limits on the spectral energy density (erg cm ^ s ^) as a function 
of gamma-ray energy for the four dSphs. 
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Fig. 2. — Exclusion regions in the (M^, {crv)) parameter space based on the results of the 
observations . It is computed according to eq. H] using a composite neutralino spectrum (see 



Wood et al\ (120081 )) and the values of J from Tabled] Black asterisks represent points from 



MSSM models that fall within ±3 standard deviations of the relic density measured in the 
3 year WMAP data set JSpergel et ndboOTh . 



